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ABSTRACT 


The  ratio  of  geostrophic  drag  coefficients  C/Cq  under 
various  stability  conditions  is  obtained  as  a  function  of 
iiichardson  number,  A  drag  coefficient  Cf  is  defined  in 
relation  to  the  mechanical  mixing  length,  and  its  ratio  to 
the  actual  drag  coefficient  C  is  related  to  the  Richardson 
number.  By  using  observational  recent  data  on  the  normal¬ 
ized  logarithmic  wind  shear  of  Monin-Obukhov  one  can  obtain 
Cf/Ca. 

The  writer  is  deeply  Indebted  to  Dr,  P.  L.  Martin 
(Professor  of  Meteorology)  for  his  suggestion  of  the  topic 
and  his  continued  help  throughout  the  investigation  and 
during  the  preparation  of  this  paper. 
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1.  Introduction. 


One  of  the  important  problems  in  meteorolo£y  is  to 
establish  relationships  betv/een  the  turbulent  characteris¬ 
tics  of  the  lower  atmosphere  and  the  large-scale  synoptic 
parameters.  A  significant  turbulent  characteristic  is  the 
ground  drag,  which  enables  one  to  specify  the  v;ind  profile, 
the  eddy  diffusivity  and  the  local  enerpj  dissipation  in 
the  surface  layer,  by  knowing  the  topographical  characteris¬ 
tics  i.e.,  the  roughness  parameter  of  a  particular  place, 
Lettau  [7]  had  made  an  attempt  to  relate  the  drag 
coefficients  with  the  static  stability  represented  by  a 
dimensionless  parameter  called  Richardson  number.  His 
invest!  ations  were  based  on  1S53  O'Neill,  Nebraska,  data. 

It  is  the  main  purpose  of  this  paper  to  verify  and  extend 
his  Investigations  using  1956  C 'Neill  data.  It  is  im.portant 
to  note  that  this  to  rk  deals  with  a  relatively  smooth 
surface,  tljerefore  the  "zero-plane  displacement"  which  is 
sometimes  used  in  the  adiabatic-wind  profile  has  been 
neglected,  a  procedure  which  greatly  simplifies  the  com¬ 
putations. 
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2.  The  Nature  of  the  Data. 


The  data  used  in  this  study  were  obtained  by  the 
Project  Prairie  Grass  conducted  near  O'Neill,  Nebraska 

during  the  summer  1956.  Since  this  study  rakes  use  of  some 
of  the  published  data,  i.e.,  v/ind  speed,  temperature  and 
surface  maps,  it  is  necessary  to  describe  some  aspects  of 
the  instrumentation.  The  instrumentation  employed  was  that 
of  the  mobile  rnicrometeorolof ical  station  of  the  Texas  A&M 
proup.  The  station  has  a  slender  aluminum  mast  supporting 
six  anemometers  at  heights  of  8,  4,  2,  1,  0,5  and  0,25  meter 
A  similar  mast  supported  seven  temperature-measuring,  radia- 
tion-shielded  copper-constantan  themocouple  junctions  at 
heights  8,  4,  2,  1,  0.5,  0,25  and  0.125  raeters.  In  addition 
the  published  report  of  Project  Prairie  Grass  includes 
sectional  sea  level  pressure  maps  at  times  nearlj’^  synoptic 
with  the  micrometeorological  data.  These  maps  revealed  some 
small-scale  features  of  the  circulation  near  O’Neill  which 
did  not  appear  on  the  iarger-scale  facsimile  maps.  The  micro 
meteorological  data  periods  v;ere  selected  to  be  nearly  simul 
taneous  with  the  available  maps.  At  times,  linear  interpola 
tion  between  the  available  charts  v/as  used  for  computations 
of  geostrophic  wind  speeds, 

Regarding  the  errors  of  measui'ernent ,  two' types  of  tem¬ 
perature  error  occur:  i.e.,  calibration  and  radiation  error. 
An  estimate  for  calibration  error  was  0.05OC.  The  estimated 
radiation  error  was  O.l^^C  with  the  understanding  tl  at,  in 
daytime  with  clear  skies,  and  low  v/indspeed,  all  measured 
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air  temperatures  overestimate  the  true  temperature.  At  night 
''/Vith  clear  skies  the  measured  temperatures  v/ould  he  lower 
than  the  ambient  temperature.  For  windspeed,  measurements 
were  dependent  upon  different  techniques  of  calibrations 
according  to  the  windforce.  All  measurements  used  in  this- 
paper  were  twenty-minute  means.  The  reader  is  referred  to 
[l]  for  more  details  concerning  the  instruinentation. 


3.  Background  Theory. 

Extensive  use  will  be  made  of  Prandtl’s  turbulence 
model,  because  of  its  relative  simplicity.  This  model  is 
described  in  many  standard  texts.  According  to  this  theory, 
the  followin  major  results  emerge: 


a  =1'^ 

dz 


X'  =  cff(ll)' 


(1) 


(2) 


where  U,'  ,  fuh' ,  are  th  e  turbulent  velocity  fluctuations  in 
the  direction  of  the  mean  wind  and  in  the  vertical,  respec¬ 
tively.  ^  is  the  eddy  mixin.'"  length,  while  ^  is  the  root- 
mean-souare  mixing  len>  th.  The  eddy  stress  X'  is  then 
considered  to  he  very  nearly  a  constant  Xq  within  a  layer 
called  the  surface  layer. 

In  the  case  of  the  fully-rough,  neutral  surface  ha^^er, 
Rossby  has  shov;n  thatXski,  so  that  (2)  may  he  integrated 
to  give 

M-O,  =  yL.  4  (3) 


where  VTa/f  is  the  friction  velocity  in  a  neutral 

surface  layer,  is  the  roi^p'hness  paraneter  and  k=0,38 


is  the  Von  Karman  constant. 

In  the  non-neutral  surface  layer,  numerous  theoretical 
profiles  have  been  derived.  However  Deacon  has  proposed  the 
semi-empirical  relationship 


^  -  Hi  f  ^ 

dz  kz^\2.oJ _ 

where  the  friction  velocity  u.*.  Is  constant  in  the 


(4) 


non-neutral  surface  layer.  The  parameter  6  is  a  decreasing 
function  of  Richardson  number.  Actually  Davidson  and 
barad  [43  have  shown  that  [0-l\  increases  somewhat  with 
height.  However  (4)  still  affords  a  useful  tool  as  an 
overall  or  bulk-relationship.  Two  useful  results  may  be 
derived  from  (4).  These  are; 

(i)  the  integrated  wind  profile; 

1-0  T 


U.  r: 


U., 


[ft)- 


(5) 

kO-0)  “  '  ' 

(li)  the  value  JU.^  by  the  llmitin  process  indicated 
in  eouation  (7)  below 


a. 


 Sm/S  ^3 


(3) 


so  that 


(7) 


Actually  the  lowest  level  at  which  Avind  data  were  available 
was  Z  =  25  cm  and  Zq  is  the  rouphness  ja  rameter  determined 
under  the  time-nearest  neutral  conditions,  may  therefore 
be  approximated  by  the  finite-difference  expression 
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— 


0. 38  AA-js" 


(6) 


jtw45-J^Z, 

Following  a  procedure  similar  to  that  employed  by  David 
son  and  Farad  [4]  in  neutral  cases,  it  may  be  shown  that 

where 


- 


*611171^ 


(9) 


i  e) 


(10) 


and 


AU.  =  (11) 

for  mean-wind  measurements  at  25,  50,  100,  200,  400,  800  cms 


4.  Geostrophic  Drag  Coefficients, 

Lettau  W  ’as  defined  the  reostropbic  drag  coefficient 
as  » 

C  =  a.  /v,  (12) 

For  neutral  stability,  the  corresponding  drag  coefficient 
is  defined  as 

Neutral  wind  profiles  at  O'Neill  may  be  obtained  by  finding 

those  cases  for  which  the  potential  temperature  profiles 

we  characterized  by30/5isO  in  the  layer  25  cm  -  200  cm. 

Altogether  15  such  profiles  were  found.  Assuming  a  lo'^arith 

mic  wind  profile  can  be  fitted  to  the  wind  data,  of 

equation  (5)  may  be  obtained  as 

Vau.  ^  •  “ 


9^2 


(14) 


By  analogy  with  flow  In  circular  conduits,  Lettau 
has  suggested  a  drag-coefficient  Ca  for  neutral  cases  of 
th  e  f o  rm 


0.104 


(15) 


”  jlo3(C«.Ro)-a.2^4 

where  ilo  may  be  called  the  "surface  Rossby  number"  defined 

as  the  non-dimensional  expression 


•Ro=V3/fz,  (IS) 

The  constants  of  eouation  (15)  were  obtained  by  Lettau  by 

least-square  methods  applied  to  a  wide  variety  of  neutral- 
wind  profiles  for  different  sites., 

Lettau  goes  on  to  give  values  of  C/Cg  for  different 

values  Richardson  number  at  100  cm,  for  the  1952 

O’Neill  data.  He  obtains  a  curve  as  shov/n  in  fig,  1, 


Ratio  of  geostrophic  drag  coefficient  for  non-adiabati c  and 
adiabatic  conditions  as  a  function  of  liichardson  number  at 
Z =  100  cm,  at  O’Neill,  Nebraska  1953  (After  Lettau  CvJ). 

It  must  be  pointed  out  tbat  here  represents  the  geostro¬ 
phic  drag  coefficient  for  any  surface  layer  of  drag  coeffi¬ 
cient  C,  if  this  layer  were  suddenly  to  be  converted  into 


a  neutral  layer  without  changing  Vg/fzQ, 

A  primary  objective  of  this  paper  was  to  test  eouation 
(15)  as  a  means  of  determining  Cg  based  on  the  1956  O'Heill 
data.  Also  application  of  non-neutral  profile  data  to  eoua- 
tions  (7)  and  (12)  afforded  values  of  C  under  identical 
roughness  conditions.  The  comparison  of  Lettau's  1953  graph 
of  C/Cq  versus  l^iioo  with  that  obtained  by  the  writer  at 
ilil40  (to  be  more  exact,  the  Richardson  number  at  141. 4cm) 
based  on  the  1956  O'Neill  data  is  shown  in  fig,  3. 

If  one  knows  a  value  of  Ro  =  Vf/fZo»  it  is  possible  to 
solve  equation  (15)  for  Cg  by  numerical  methods.  It  is 
shown  in  the  Appendix  that  a  numerical  solution  Cg  of  eaua= 
tlon  (15)  may  be  obtained  iteratively  from  the  eouation 

V,-(a.)024  +  y- +  6.5-927  _  ) 

V  s  'it  -  — - - ^ - - - — -  (17) 

where  yi=  i-th  iterant  to  0.2395/Cgj,  The  first  estimate 


7 


The  chief  inaccuracy  in  this  approach  appears  to  be  uncer¬ 
tainty  in  the  numerical  constants  appearing  in  equation ( 15) , 
especially  at  low  geostrophic  v;lnd  speedso  Another  source 
of  error  is  the  necessity  in  some  cases,  of  interpolating 
Vg  between  two  maps  6  or  12  hours  apart  in  order  to  get  the 
estimated  value  of  Vg  at  the  time  of  micrometeorological 
data-run. 

A  second  approach  to  determine  the  relationship  betv/een 
C  and  Cq  was  used.  Based  on  15  neutral  cases  at  O'Neill, 
Nebraska,  values  of  v/ere  computed  using  equations  (13) 
and  (14).  Values  of  Cg  were  then  plotted  versus  log  Ro, 
and  the  resulting  best-fitting  curve  compared  to  Lettau's, 
which  is  replotted  on  the  same  graph.  The  comparison  will 
be  discussed  in  later  section. 

The  liichardson  number  has  been  used  very  extensively 
as  a  measure  of  stability  in  the  micro-meteorology.  Its 
definition  is 


j  bG/d£ 


(18) 


0  (  5SL/diy 

This  number  is  proportional  to  the  ratio  of  work  done  against 
static  stability  to  the  work  done  by  eddy  stress. 

For  instruments  arranged  at  eoual  logarithmic  spacing, 
as  in  the  1953  and  1956  O'Neill  data,  a  convenient  finite- 
difference  form  of  Ri  at  Z=:Z]^^2  been  proposed  by 

Lettau  [sj 

i-w/  A  ch  9.  ^  ( o.  ~  §.)  ]Lt^  0. 

(19) 


p./_  S  Zut  ^  2. 


This  value  is  understood  to  be  applicable  at  the  geometric 
of  Zp  and  Z2,  namely  0  2  and  0-2  in  the 
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mean 


1956  data  are  the  20-ininute  mean  potential  temperature  and 
v/indspeed  at  level  2,  and  an  analogous  meaning  applies  for 
symbols  with  subscript  one. 

Finally,  Monin  and  Obukhov  have  introduced  a  use¬ 
ful  non-dimensional  parameter  S  defined  by 

r  _  kz  ad/bZ  (2Q 

Note  that  the  nximerator  of  S  has  the  form 


ki  SCL/9z  r:  if 


(21) 


which,  according"  to  the  Prandtl  mode],  may  be  regarded  as 
representing  here  defined  represents  that 

part  of  UL,,  due  to  turbulent  fluctuations  associated  with 
the  mechanical  turbulent  raixin”  lenpth  ifakl.  The  expression 
for  3  may  then  be  written 


C.  —  ^»f  / ^3 

^  ■  a,/V3  -  c 


The  drap  coefficient  defined  by  Cf  =  U.n^/Vg  is  not  one 
that  can  be  measured  directly,  but  is  one  whose  value  can  be 
Inferred  by  empirical  information  repardinp  S,  This  is 
dealt  with  in  more  detail  in  section  8,  This  interpretation 
of  as  kz  3u./2>Z  appears  to  have  been  piven  first  by 

Businper  [s}  , 


5,  Comparison  of  C/Cq  for  1952  and  1956  data. 

In  computing  Cg,  Lettau's  empirical  formula,  ecuation 
(17)  was  first  used  with  the  1956  data.  It  was  necessary 
to  compute  Zq,  usinp  ecuation  (9)  from  the  time-nearest 
neutral  wind  prol'iles  and  Vp  from  interpolated  surface 
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charts.  In  ccxputing  U.-*and  C,  eouation  (8)  and  (12)  were 
used.  The  value  of  Rii40  is  also  included  using  equation 
(19).  The  results  for  all  data-days  employed  are  given  in 
table  1,  and  the  groupings  ol  ilii40  classes  are  given  in 
table  2.  The  graph  of  C/Ca  versus  Rii40  is  given  in  fig.  3, 
with  Lettau's  curve  hi  superimposed.  The  agreement  with 
Lettau's  results  is  not  as  good  as  might  have  been  expected. 
For  instance  the  value  C/Cg^~  1,75  was  found  at  Ri]_4Q-0 
contrasted  with  Lettau's  value  C/Cq-I.OO  at  Rl]_oO"®‘>  More¬ 
over  the  maximum  value  of  C/Ca  found  in  this  paper  exceeds 
Lettau's  maximum  by  a  comparable  factor  but  with  RI140  ■= 
-0.055,  compared  with  Hli00=  -0.007,  according  to  Lettau, 
Possible  sources  of  error  lie  in  the  values  of  Vg, 
which  affect  the  values  of  C  and  Cg  given  by  equations  (12) 
and  (15)  respectively,  as  well  as  in  the  value  of  Zq,  The 
values  of  Zq  determined  here  ranged  from  0.35cm  to  3,1  cm, 
as  contrasted  with  values  quoted  by  Blackadar  et  al  [2]  of 
Zo=  0.6  cm  for  profiles  in  approximately  the  same  data- 
period. 
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TABLE  I  -  Values  of  C/Cq 
Date  CST  (oms)  U.^^^m/sec) 

10  Jul  1956 


1305 

1.10 

24.67 

1405 

1.10 

35.85 

1505 

1.10 

31.24 

1530 

1.10 

30.99 

1605 

1.10 

35.75 

1705 

1.10 

35.00 

1905 

1.10 

44.00 

2105 

1.10 

21.03 

Jul 

1956 

0930 

3.14 

49.26 

1105 

3.14 

42.85 

1205 

3.14 

46.69 

1305 

3.14 

52.92 

1505 

3.14 

52.55 

1705 

3.14 

59.  69 

1805 

3.14 

67.38 

2130 

3.14 

31.03 

2305 

3.14 

31.86 

i  Jul 

1956 

0205 

1.49 

2.15 

0405 

1.49 

12.66 

0805 

1.49 

50.49 

0905 

1.49 

52.24 

1230 

0.70 

58.87 

1405 

0.70 

65.24 

1605 

0.70 

62.80 

2530 

0.70 

48.72 

'  Aup 

1956 

1805 

0.35 

29.53 

2105 

0.35 

11.65 

2205 

0.35 

25.00 

equations  (3),  (12),  (15). 


^ts) 

C 

Ca 

C/Ca  ^140 

60 

0.036 

0.034 

1,08  -0,46 

,80 

0 . 056 

0.034 

1,65  -0,17 

00 

0.052 

0.034 

1.52  -0,28 

64 

0.040 

0.033 

1.20  -0,25 

60 

0.053 

0,034 

1,54  -0,20 

80 

0.065 

0,035 

1,87  -0,05 

52 

0.084 

0,035 

2,41  -0,02 

CO 

0.035 

0,034 

1,02  0,06 

44 

0.051 

0,037 

1,38 

-0,15 

00 

0.045 

0,037 

1,22 

-0,21 

00 

0,052 

0,037 

1,39 

-0,16 

00 

0.062 

0,038 

1,66 

-0,09 

20 

0,065 

0,038 

1,72 

-0.14 

00 

0.075 

0,038 

1.97 

-0,06 

80 

0.085 

0,038 

2.25 

-0,02 

52 

0,059 

0.040 

1,47 

0,04 

52 

0.085 

0.042 

2.01 

0,03 

80 

0,009 

0.040 

0,22 

0,10 

60 

0,026 

0.037 

0,72 

0,05 

00 

0.042 

0,034 

1,22 

-0,10 

00 

0,040 

0.032 

1,24 

-0,01 

00 

0,038 

0,029 

1,27 

-0,02 

00 

0.028 

0,028 

1,02 

-0,02 

00 

0,023 

0,027 

0,85 

-0,01 

80 

0,022 

0,028 

0,78 

0,00 

60 

0.056 

0.030 

1.84 

0,00 

80 

0,022 

0,030 

0,71 

0,07 

00 

0,045 

0,030 

1,50 

0.01 

by 

Vf 

15 

12 

12 

15 

11 

10 

10 

12 

19 

19 

18 

17 

16 

16 

15 

9 

7 

4 

9 

24 

26 

32 

46 

54 

44 

10 

10 

11 
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TABL^j  1  -  (Cont’d) 


Date  CST  Zq  (cms)  a^^pm/sec) 
7  Aup  1956 


0405 

3.06 

57.72 

0805 

3.06 

63.14 

0905 

3.06 

109.10 

1105 

3.06 

69.65 

1530 

0.  62 

59.79 

1605 

0.62 

64.83 

1805 

0.62 

36.02 

2205 

0.62 

23.87 

2205 

0.62 

9.47 

27 

Auf  1956 

1230 

0.55 

42.96 

1505 

0.55 

38.09 

1705 

0.55 

40.08 

1930 

0.55 

1,79 

2205 

0.55 

2.98 

29 

Aug  1956 

2135 

1.53 

24.09 

Vg  (kts) 

C 

Ca 

C/Cg  x81i40 

11.88 

0.097 

0,039 

2,47 

0,01 

12,88 

0,098 

0,039 

2.52  -0,04 

13.08 

0.167 

0,039 

4,30  -0,02 

13.43 

0.103 

0.039 

2.67  -0,08 

15.84 

0,075 

0,031 

2,44  -0,04 

15.00 

0,086 

0,031 

2,77  -0,01 

15.76 

0,046 

0.031 

1,47  -0,01 

14.00 

0.034 

0.031 

1,09 

0,05 

12.40 

0,015 

0.032 

0,48 

0,13 

26.48 

0,032 

0.029 

1,12 

-0,07 

24.80 

0.031 

0.029 

1,06 

-0,02 

25.20 

0.032 

0,029 

1,10 

-0,02 

15.92 

0,002 

0,031 

0.07 

0,15 

26.  60 

0.002 

0,029 

0.07 

0.04 

7.20 

0.067 

CO 

o 

o 

o 

1,75 

0,03 

TABLE  2  -  Class  means  of  Hi]_40  vs  C/Cg 


^^140 

-0.43 

-0.37 

-0.31 

-0.25 

-0.19 

-0.12 

-0.06 

0,00 

0.06 

0,13 

C/Ca 

1.08 

- 

1,52 

1,20 

1,45 

1,44 

2,04 

1,72 

1,11 

0,26 

12 


Pi£.  3 

i^atio  of  geostrophic  drag  coefficient  for  non-adiabatic  and 
adiabatic  condi-tions  as  a  function  of  Richardson  nvimber 
at  Z  = 140  cm,  at  O’Neill,  Nebraska  1956. 
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6.  Revised  formula  for  Cg 

Hitherto  Gg  has  been  computed  by  usinc’  Lettau’s  empirical 
function  for  Cg,  eouatlon  (15).  Since  Lettau  himself  realized 
that  the  namerlcal  constants  in  (15)  were  not  very  exact  [?]  , 
an  attempt  is  made  here  to  obtain  an  improved  relationship 
between  Cg  and  Ro  by  uslnr  only  1956  O'Neill  data.  The  first 
step  again  consists  in  computing  Zq  for  the  time-nearest 
reutral  wind  profile  according  to  eouation  (9)o  As  beforCj, 
Uifcot  results  from  equation  (14)  with  k=0o38,  so  that 


o.o6i  AU.O.  ' 

Eouation  (13)  is  then  applied  to  determine  Cg.  The  results 
of  the  computations  are  shown  in  table  3. 

An  empirical  function  for  Cg  was  then  obtained  by  the 
following'  procedure: 

Let  log  Ro-x  and  Cgts  y  and  suppose  the  desired  function  has 
the  form 

^  =  be  (24) 

where  y  means  the  desired  function,  while  b  and  eg  are  cop- 
stants  to  be  determined.  Taking  the  log  of  (24)  one  obtains 


^25) 

To  simplify  the  notation,  let  !L>.y  and  ivb  scw,  and  there= 
fore  the  desired  function  becomes 

Y  =  O-  -V  c)£-  (26) 

Using  the  least-squares  method,  the  normal  ecuatlons  are 

Y-ol-ck.  =  o 

^  nv  ( on  \ 

-a2«,-  -cIxUo 

...  is! 

Substituting  the  values  of  observations  (table  3)  in  the 
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normal  equations  one  obtains 


3.51  -  a  -  7.21c  =0 


38L90-10S.28a  »  785.180  =.0 


from  which,  it  follows  that 


a=-0.85  or  b^O.427 


c=  -0.359 


Substitutinf  these  constants  In  (24),  and  further  by  chanf= 
inf  y  to  Cq^,  x  to  lof  Ro,  the  expression  for  c^n  be 
obtained 


(28) 


The  fraph  of  eouation  (28)  is  shown  In  fip,  4.  compared  with 
Lettau’s  praph.  From  fif.  4  it  is  evident  that  the  revised 
formvila  for  Cq  tends  to  underestimate  the  draf  at  larfer 
values  of  Rossby  number,  and  overestimate  it  at  smaller  Ro 
(lof  Ro<7.1).  It  should  be  noted,  however,  that  there  was 
considerable  scatter  of  the  observations  relative  to  the  line 
of  best  fit. 

If  one  assumes  that  ecuation  (28)  is  a  universal  formula 
for  eostrophic  draf  coefficient  in  niddle  and  higher  latitudes 
under  neutral  conditions,  one  can  then  write 


(29) 


with  Cq  flven  by  (28).  Hence,  adiabatic  wind  speed  at  level 
Z  is  obtainable  as  a  function  of  Vf,  Zq,  Z  and  latitude. 
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TABLE  3«  ”  Results  of  computations  for 


Adiabatic  hours 

2q  cm 

UL*^6ra/  sec) 

Vg  (jets) 

yi  =  Ca 

Xi=.log  Ro 

li  =  lr>l 

Yi  xi 

10  July  3956 

1905  CST 

1.10 

36.66 

11.00 

0.067 

6.707 

-2.703 

-18.1 

11.89 

23  July  1956 

1905  CST 

3. lit 

31.10 

11.00 

O.OI5 

6. 355 

-3.101 

-19.1 

10.15 

25  July  1956 

0605  CST 

1805  CST 

l.ii9 

0.70 

29.59 

57.31 

16.80 

61.00 

0.035 

0.016 

6.759 

7.670 

-3.352 

-1.017 

-22.6 

-30.7 

15.70 

58.83 

6  Aug  1956 

1905 

0.35 

21.79 

10.60 

O.OI6 

7.192 

-3.079 

-22.1 

51.70 

7  Aug  1956 
n6o5  CST 

1905  CST 

3.06 

0.62 

36.12 

32.91 

12.20 

ll.OO 

0.056 

o.oll 

6.313 

7.08I 

-2.817 

-3.121 

-17.9 

-22.1 

39.62 

50.13 

27  Aug  1956 

1705  CST 

0.55 

10.69 

25.20 

0.032 

7.368 

-3.112 

-25.1 

51.32 

29  Aug  1956 

1930  CST 

1.53 

36.72 

32,00 

0,032 

8.025 

-3.6-17 

-30,6 

61.18 

11  July  1956 

06o5  CST 

O.li 

11.61 

ll.OO 

0,020 

7.618 

-3.912 

-30.6 

61.62 

?h  July  1956 

o6o5  csr 

1905  CST 

0.ii9 

0.61 

27.57 

23.79 

20.00 

21.20 

0.028 

0.022 

7. 131 
7.121 

-3.576 

-3.817 

-26.5 

-27.1 

55.20 

50.69 

26  July  1956 

1605  CST 

0.61 

37.52 

lo.oo 

0.019 

7.398 

-3.963 

-29.3 

51.76 

6  Aug  1956 

1905  CST 

0.3? 

32.61 

32.00 

0.020 

7.615 

-3.912 

-29.8 

58.52 

26  Aug  1956 

1P05  CST 

Sums 

Average 

0.71 

27.57 

32.00 

0.017 

7.360 

108.279 

7,210 

-1.075 

-52.737 

-3.51 

-30,0 

-381.9 

51.17 

785.16 

Pig.  4 

Geostrophic  drag  coefficient  under  neutral  conditions  as 
a  function  of  surface  Rossby  number. 
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7. 


A  revised  formula  for  C/Cg 

The  revised  formula  for  Cg  will  be  used  In  this  section 
to  compute  C/Cq,  still  usinf  the  previously  tabulated  values 
of  peostrophic  v/inds  and  C  (table  3)o  The  results  are  shown 
in  table  4  and  tie  resultinr  values  cf  C/Cg  then  plotted 
versus  Rii4o  as  shown  in  fif.  5.  The  procedure  of  fittinp 
the  curve  Is  shov;n  below. 

After  plotting  the  values  of  C/Cg,  a  further  restriction 
is  placed  on  the  empirical  function  C/CaJ  namely  that  C/Cgsl 
at  Ri=0.  ivoreover,  since  the  shapes  of  the  curves  for  C/Cjj 
in  unstable  and  stable  conditions  are  different  (according 
to  the  scatter  cf  points),  it  is  assumed  that  in  stable  condi¬ 
tions  the  desired  function  has  the  form  y-i  where  y^C/Cg 

and  X  r  *iip4o*  unstable  cor.diticn,  y  is  assured  to  be  of 

k  X. 

the  form  \  +  *  where  yzC/Cg  and  x-  “Rii40i>  a,  kp, 

k2,  are  constants  to  be  deterr  ined, 

(i)  For  stable  conditions  (table  5),  with  Rip4o=x  and 
C/Cgs  y,  if  the  assumed  function  has  the  form 


(30) 


then 


!U  y  ~ 


and  if  =  Y 


V  -  ctic 


Bv  the  method  of  Icast-s cuares ,  th-  normal  eouation  will  be 
«/  ' 


(31) 


Solvinf  for  a. 


^ j  S  Q 
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Sutstitutlnc  this  constant  in  (30)  and  chan  inf  y  to  C/Cg 
X  to  Rii40>  the  expression  for  C/c^  is  obtained  as 


C/C<x=  ^‘‘^9 


(32) 
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TABLE  Ij.  -  Computations  of  C/C^  using  equation  (28)  for 


Date  CST 

log  Ro 

r- 

Ca 

C 

10  July  19^6 

1305 

6.799 

-2.51 

0.081 

0.03h8 

0.0363 

i.oh5 

-0.h6 

1U05 

6.773 

-2.50 

0.082 

0.0350 

0.0560 

1.600 

—0 . 17 

1505 

6.7li5 

-2.1(9 

0.083 

0.0360 

0.0521 

i.hh5 

-0.28 

1530 

6.860 

-2.53 

0,080 

0.03h2 

0.0396 

1.160 

-0,25 

l6o5 

6.730 

-2.1,8 

O.OBh 

0.0359 

0.0530 

1.1(75 

-0.20 

1705 

6.699 

-2.1i7 

0.085 

0.0363 

0.06h8 

1.785 

-0.05 

l8o5 

6.687 

-2.U6 

0.085 

0.0363 

0.0838 

2.320 

-0.02 

2105 

6.715 

-2.1(9 

0.083 

0.0360 

0.0350 

0.975 

0.06 

23  July  1556 

0930 

6.500 

-2.  ho 

0.091 

0,0388 

0.0507 

1.305 

-0.15 

1105 

6.ii87 

-2.39 

0.092 

0.0393 

0.01(51 

l.lh5 

-0.21 

1205 

6.li6ii 

-2.38 

0.093 

0.0397 

0.0519 

1.305 

-0.16 

1305 

6.h39 

-2.37 

0.093 

0.0397 

0,0623 

1.570 

-0.09 

l6o5 

6.ia8 

-2.36 

0.09h 

O.OhOl 

0.06h9 

1.615 

-O.lh 

1705 

6.1(13 

-2.36 

0.09h 

O.OhOl 

0.07h6 

1.860 

-0.06 

1805 

6.I107 

-2.36 

0.09h 

O.OhOl 

0.0853 

2.150 

-0.02 

2130 

6.230 

-2.30 

0.100 

O.Oh27 

0,0591 

1,3B‘^ 

O.Oh 

2305 

s.oeh 

-2.22 

0.109 

O.Oh65 

0.0850 

1.830 

0.03 

25  July  1956 

0205 

6.215 

-2.29 

0.101 

0.0h31 

0.0090 

0.209 

0.10 

0U05 

6.516 

-2.  ho 

0.091 

0.0386 

0.026h 

0.680 

o.o5 

o8o5 

6.721 

-2.h7 

0.085 

0.0363 

0.0h21 

1.160 

-0.10 

0905 

6.9U8 

-2.56 

0.077 

0.0329 

0.0h02 

1.220 

-0.01 

1230 

7.366 

-2.72 

0.066 

0.0282 

0.0368 

1.305 

-0,02 

Ili05 

7.523 

-2.77 

0,063 

0.0269 

0.028h 

1.055 

-0.02 

1605 

7.593 

-2.80 

0.061 

0.0260 

0.0233 

0.895 

-O.Gl 

2330 

7.512 

-2.77 

0.063 

0.0269 

0.0218 

0.810 

0.00 

6  Aug  1956 

1805 

7.188 

-2.65 

0.071 

0.0303 

0.0557 

1.835 

0.00 

2105 

7.196 

-2.66 

0.070 

0.0299 

0,0216 

0.722 

0.07 

2205 

7.20l( 

-2,66 

0.070 

0.0299 

0.0h5h 

1.515 

0.01 

7  Aug  1956 

0ij05 

6.295 

-2.32 

0.098 

0.0hl8 

0.0970 

2.320 

0.01 

0805 

6.330 

-2.33 

0.097 

0.0hl5 

0.0979 

2.360 

-O.Oh 

0905 

6.337 

-2.3h 

0.096 

o.ohio 

0,1666 

h.060 

-0.02 

1105 

6.350 

-2.3h 

0.096 

O.OhlO 

0.1032 

2.510 

-0.08 

1330 

7.111 

-2.62 

0.073 

0.0312 

0.0755 

2.h20 

-O.Oh 

1605 

7.087 

-2.61 

0.07h 

0,0316 

o,o86h 

2.730 

-0.01 

i8o5 

7.108 

-2.65 

0.071 

0.0306 

O.Oh37 

l.^hO 

-0.01 

2005 

7.057 

-2.60 

0.07h 

0.0316 

0.03hl 

1.069 

o.o5 

2205 

7.000 

-2.56 

0.07h 

0.0329 

0.0153 

o.h65 

0.13 
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TABLE  it.  (Cont'd) 


»»s 


Date  CST 

log  Ro 

c  C/Ca 

^iliiO 

27  Aug  1956 

1230 

7.390 

-2.73 

0.065 

0.0278 

0.0325  1.170 

-0.07 

1505 

7.361 

-2.72 

0.066 

0.0282 

0.0307  1.090 

-0.02 

1705 

7.368 

-2.72 

0.066 

0.0282 

0.0318  1.125 

-0.02 

1930 

7.168 

-2.6U 

0.071 

0.0303 

0.0023  0.076 

0.15 

2205 

7.392 

-2.72 

0.066 

0.0282 

0.0022  0.078 

o.oU 

29  Aug  1956 

2135 

6.377 

-2.35 

0.095 

0.0ij06 

0.0669  1.6ii0 

0.03 

TABLE  5. 

C/C^i  in 

stable  conditions 

• 

L 

C/Ca^y; 

M  »! 

1 

0.057 

0.975 

-0.025 

-0.001 

O.OC3 

2 

0.036 

1.385 

0.326 

0.012 

0.001 

3. 

0.035 

1.830 

0.6oit 

0.021 

0.001 

h. 

0.099 

0.209 

-1.565 

-0.155 

0.010 

5. 

0.051 

0.680 

-0.386 

-0.020 

0.003 

6. 

0.067 

0.722 

-0.326 

-0.022 

o.ooil 

7. 

O.Olit 

1.515 

o.i4i5 

0.006 

o.oco 

8. 

0.008 

2.320 

0.81t2 

0.007 

0.000 

9. 

o.o5ij 

1.069 

0.068 

O.OOit 

0.003 

10. 

0.127 

O.i.65 

-0.766 

-0.010 

0.016 

11. 

0.1'^3 

0.076 

-0.27i| 

-0.0i)2 

0.023 

12. 

O.OitO 

0.078 

-0.2it8 

-0.010 

0.002 

13. 

0.031 

1 .6ii0 

o.it95 

0.015 

0.001 

UlTiS 

-0.283 

0.066 
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(ii)  For  unstable  conditions  (table  6),  the  scatter 
of  points  of  C/Cq  versus  Rii40  suggests  a  curve  of  the 
form 

izz) 

where  -  Rii4c  =  x  and  C/C^-z.  y.  Again  C/Cq  has  been  forced 
to  approach  unity  as  Ri  approaches  zero.  If  we  let  y-l-'f 
and  then  take  natural  logs 

^  If  s  £vt  k,  +  X,  +  k,«.  (34) 

SUt  =  ♦  U 

Let  Ll-.\  and  JL/  k,  S  ct  ,  and  substitute  in  (34) 

Y  =  A.  4-  kjX, 

The  normal  eouations  are 


Y  -  a  -  o 

rw 

2Yx,.- 

(W  ^ 

iZl 

Solving  for  and  k2 


kir  30.364  ,  kg  =  -13.324 

Substituting  these  constants  in  (33)  and  replacing  x  by 
_Ril40,  y  by  C/Cg,  gives 

C/C^  =  «  -  ^’.40  Sitp  [l3.32A‘Rt,4o] 

Note  that  C/Cg  of  equation  (36)  may  be  maxiniized  for  Rii40 
-  ‘  -  -0,075,  for  which  (C/Ca)j„gx  is  1,339,  This  is 

close  to  the  maximum  of  the  fitted  curve  of  fig,  5,  Lettau's 
graph  of  C/Cg  has  teen  included  in  fig.  5  and  it  may  be  noted 

that  the  computations  of  t'^io  paper  are  in  general  agreement 
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v.'ith  those  of  Lettau,  although  his  data  did  not  span  as 
a  range  of  Ri  values.  The  results  indicate  that  C/C^y  1 
stable  conditions,  whereas  C/Cg  <  1  in  stable  conditions 


large 
in  un- 
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L 

1 

2 

3 

a 

5 

6 

7 

r 

9 

10 

D1 

12 

13 

111 

15 

l6 

17 

18 

19 

20 

21 

22 

23 

2h 

25 

26 

Su-TTI 


TABLE  6  -  C/Ca  in  Unstable  Conditions 


Yu; 


0.ii63 

0.165 
0.280 
0.251 
0.199 
0.053 
0.017 
0.11i7 
0.210 
G.]^7 
0.091 
O.liil 
0 .058 
0.022 
0,103 
0.011 
0.017 
0.019 
0.01^2 
0.077 
0.0h3 
0.011 
C.OlO 
0.075 

0.02ij 

0.02U 

2.710 


Vi  =c/c^ 

1.0U5 

1.600 

l.ijli5 

l.l6o 

1. klS 
1.785 
2.320 

1.305 
1.1^5 

1.305 
1,570 
1.6i5 
1.860 
2.150 
l.l6o 
1.220 

1.305 
1.055 
2.360 
2.510 

2. U20 
2.730 
1.5U0 
1.170 
1.090 
1.125 


0.0h5 

0.600 

oMS 

0.160 

0.ij75 

0.785 

1.320 

0.305 

0.1U5 

0.305 

0.570 

0.615 

0.860 

1.150 

0.160 

0.220 

0.305 

0.055 

1.360 

1.510 

1.U20 

1.730 

o.5bo 

0.170 

0.090 

0.125 


0.097 

3.630 

1.590 

0.638 

2.390 

IU.700 

77.600 
2.080 
0.690 

1.9i40 

62.600 
U.360 

111.  800 

52.3  DO 
1.550 
20.000 
17.950 
2.890 
32.1^00 

19.600 
33.000 

157.300 

5I4.OOO 

22.700 

3.750 

5.200 


-2.333 

1.289 

O.hSh 

-Oj'h9 

0.871 

2.688 

I4.352 

0.737 

-0.371 

0.663 
a. 137 

l.h72 

2.695 

3.957 

0.ij38 

2.996 

2.888 

1.061 

3.1^78 

2.976 

3.500 

5.058 

3.989 

3.122 

1.322 

1.6^9 

52,6ii9 


O.21I1 

0.027 

0.078 

0.063 

O.OUO 

0.003 

0 

0.022 

O.OUi 

0.025 

0.008 

0.020 

0.003 

0 

0.011 

0 

0 

0 

0.002 

0.006 

0.002 

0 

0 

0.006 

0.001 

0.001 

0.576 


-1.080 

0.213 

0.130 

-0.113 

0.173 

0.11^2 

0.07il 

0.108 

-0.078 

o.ioii 

0.376 

0.208 

0.156 

0.087 

o.oli5 

0.033 

0.0h9 

0.020 

0.lli6 

0.229 

0.151 

0.056 

o.oho 

0.23h 

0.032 

0.0140 

1.575 


24 


Pig.  5. 


Ratio  C/C^  as  a  function  of  Richardson  number  at  Z = 140  cm, 
at  O'Neill,  Nebraska  1956  with  given  by  equation  (28). 


8.  The  normalized — logarithmic  wind  sheai  S 

It  was  explained  in  section  4  that  if  one  knows 
S  ^  ^^/c)2  ^  jjjgy  obtain  directly  Cf/Cjor 
Ellison  [2"]  has  fiven  a  prediction  equation  for  S  in  unstable 
cases 


=  i  (37) 

where  Z  is  heipht  and  L.‘  is  the  {gradient  length  defined  as 


L‘  =  - 


a*e 


(56) 


with  defined  as 


K, 


7'-  f  (39) 

here  T  is  a  constant,  and  l<^  and  are  eddy  diffusivity 

coefficients  for  heat  and  momentum  respectively o  Dividing 
both  sides  of  (37)  by 

1  +  ^  =  5''' 

Substituting  (20)  and  (18)  in  the  expression  for  and 

solving  for  S 

5  =  .,0) 

Blackadar  et  al  [2]  suggest  the  value  ^''=18  for  urstaLle 
conditions.  However  in  stable  conditions,  they  find  that 
equation  (40)  ( ives  no  tetter  approximation  than  the  L.onin- 
Obukhov  " log +•  linear"  wind  profile  defined  by 


S  -  »  +  -21  z  c<  -consta-n^ 

"  L'  ^ 

which  is  ordinarily  valid  for  Richardson  numbers  near  cneo 
Again,  it  may  be  verified  that  (41)  is  equivalent  to  an 
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ecuation 


5-  »  4-oC'i^'  (42) 

Blackedsr  et  al  [2]  have  tabulated  values  of  t><.'  at  C’lelll 
for  various  layers.  For  the  layer  1-2  m  under  stable  con¬ 
ditions,  the  best-fitting  value  of  the  constant  was 
which  v.'ill  be  used  for  the  stable  cases  in  this  paper® 

The  computations  for  S  are  shown  in  table  7®  The  values 
of  5  are  presumed  applicable  at  the  heirht  assigned  to  Ri, 
that  is  at  Z  -  or  Z-^141,4  cm®  The  curve  of  S  is 

shown  in  fig.  6.  The  result  indicates  that  S  < 1  in  unstable 
cases,  and  S >  1  in  stable  cases.  Whether  S  in  stable  condi¬ 
tions  levels  off  with  a  further  increase  of  Hi,  is  not  cer¬ 
tain  yet,  since  the  value  chosen  for  may  not  be  valid  for 
larf  e  positive  values  of  Ri. 
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Results  of  computations  for  S 


TABLE  7.  - 


Date  CST 

^140 

Vg  (kts'j 

^(140) 

10  July  1956 

1S05 

-0.46 

13.60 

0.572 

1405 

-0.17 

12.80 

0.708 

1505 

-0.28 

12.00 

0.638 

15S0 

-0.25 

15 . 64 

0.652 

1605 

-0.20 

11.60 

0.683 

1705 

-0.05 

10.80 

0.850 

1805 

-0.02 

10.52 

0.935 

2105 

0.06 

12.00 

1.515 

23  July  1956 

0930 

-0.15 

19.44 

0.724 

1105 

-0.21 

19.00 

0.676 

1205 

-0.16 

18.00 

0.715 

1305 

-0.09 

17.00 

0  784 

1605 

-0.14 

16.20 

0.729 

1705 

-0.06 

16.00 

0.836 

1805 

-0.02 

15.80 

0.918 

2130 

0.04 

9.52 

1.326 

§305 

0.03 

7.52 

1,313 

25  July  1956 

0205 

0,10 

4,80 

1.888 

0405 

0.05 

9.60 

1.457 

0805 

-0.10 

24.00 

0.769 

0905 

-0.01 

26 . 00 

0,958 

1230 

-0.02 

32.00 

0.936 

1405 

-0.02 

46.00 

0.929 

1605 

-0.01 

54.00 

0.9f^0 

2330 

0.00 

44.80 

1.033 

G  Aut'  1956 

1605 

0,00 

10 . 60 

0.994 

2105 

0.07 

10.80 

1.605 

2205 

0.01 

11.00 

1.123 

7  Aur  1956 

0405 

0.01 

11.88 

1.076 

0805 

-0.04 

12.88 

0.868 

0905 

-0.02 

13.08 

0.939 

1105 

-0.08 

13.48 

0.804 

1330 

-0.04 

15.84 

0.867 

1605 

-0.01 

15.00 

0.958 

1805 

-0.01 

15.76 

0,960 

2005 

0.05 

14.00 

1.490 

2205 

0.13 

12.40 

2.147 
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TABLE  7.  (Cont'd) 


Date  GST 

^140 

Vr  (kts) 

^(140) 

27  Auf  1956 

1230 

-0.07 

26.48 

0.808 

1505 

-0.02 

24.80 

0.915 

1705 

-0.02 

25.20 

0,914 

1930 

0.15 

15.92 

2,375 

2205 

0.04 

26.  60 

1.362 

29  Au‘‘'  1956 

2135 

0.03 

7.20 

1,278 

29 


The  normalized  logarithmic  wind  ahear  S  as  a  function  of 
iUchardson  number  at  Z=  140  cm  based  on  1956  O'Neill, 
Nebraska  data. 
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9,  An  alternative  method  of  determination  of  '-ta¬ 
in  section  7,  C/Cg  v.’as  computed  by  eouatlons  (8),  (12), 
and  (28).  Yi/hile  there  may  be  considerable  error  in  assuming 
a  functional  form  for  of  the  type  derived  in  section  6, 
there  is  also  a  possibility  of  error  in  the  computation  of 
by  equation  (8).  Hence  an  alternative  computation  of 
was  employed,  one  which  employs  the  expressions  for  S  of  the 
section  8.  This  approach  essentially  leads  to  a  value  of 
in  terms  of  the  wind-profile  U.C5)  and  the  ir.terral-mean  value 
of  S  in  a  sub-layer  of  the  surface  layer.  Thus  for  example 

and  intefratinp  from  Z  =  Zq  to  Z  -  Z,  assuming  is  constant 

in  the  surface  layer  leads  to 


where  3  is  the  intef ral-m.ean  value  of  S  up  to  level  Z. 

As  in  the  previous  section  S  was  taken  in  accordance 

with 

3  \  unstable  cases 

U4) 

■S”  *  9  stable  cases 

Although  Priestly  has  raised  a  question  re{  ardinf  the 

the  validity  of  a  constant  with  hei;  ht  in  stable  condi«= 

tions,  it  was  assumed  that  a  sxirface  layer  at  the  depth  of 
at  least  2m  exists  for  both  stable  and  unstable  cases. 

Values  of  S  accordinp  to  eouation  (44)  were  plotted 
usinf  Hi  (Z)  for  values  of  Z-  5,  35,  70,  140,  280cr)S.  Ihe 
lowest  value  of  heipht  corresponds  to  a  level  at  the  eonetrlc 
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mean  of  the  roughness  parameter  Zq  and  Z =  25  cm,  usinp  the 
mown  values  of  and  for  this  layer.  The 

fraphs  of  S  versus -2^  Z  was  plotted  and  the  inter ral-mean 
value  S  extracted.  Fly.  7  shows  an  example  of  an  unstable 
and  stable  case  which  occurred  at  1905  and  2105  GST,  respec¬ 
tively  on  10  July  1956.  This  permits  one  to  solve  for  '-t,  by 
means  of 

_kaC5££)^  (45) 

S  (aoo/Zt,) 

The  values  of  computed  by  equation  (45)  are  displayed  in 

table  8  alonr  v/ith  the  values  of  listed  earlier  in 

table  1.  The  absolute  average  error  is  approximately  10^, 
although  the  algebraic  error  is  Z%. 

Note  tiiat  the  ratio  Cf/Cg  may  be  obtained  as  the  product 

Cc  '  c  Co. 

—  ^  (43) 

Co. 

The  value  of  S  as  a  function  of  Ri(Z)  may  be  obtained  from 
eauation  (44),  Values  of  U.^.  may  be  obtained  from  equation 
(45);  and  Cg  may  be  determined  from  equation  (28), 

The  main  purpose  of  this  last  discussion  is  to  empha¬ 
size  that  Cf  is  not  equal  to  Cg,  Listed  below  in  table  8 
are  some  computed  valu.es  of  Cf/Cg, 
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The  nomalized  logarithmic  wind  shear  S  as  a  function  of 
for  unstable  and  stable  conditions  for  10  July  1956. 


TABLE  8 


Cf/Cg  in  unstable  and  stable  cases 


Date  CST 

i^il40 

S 

sec) 

(C  f /C  a  W^sec)  from  ( 

10  July  1956 

1505 

-0.46 

0,572 

25,1 

0,605 

24,7 

1405 

-0.17 

0,900 

55,0 

1,106 

.505 

-0.28 

0,685 

51,3 

0,925 

31,2 

1605 

-0.20 

C .  900 

31,6 

1,038 

55,8 

1705 

-0.05 

0,950 

34,8 

1,515 

35,0 

1805 

-0.02 

0.97C 

43,0 

2,100 

44,0 

23  05 

0.06 

1,150 

20,5 

1,455 

21,0 

23  July  1956 

1105 

-0,21 

0,850 

55.7 

0,648 

42,9 

1205 

-0.16 

0,880 

38,7 

0,773 

46,7 

1505 

-0.09 

0,685 

44,3 

1,025 

52,9 

1605 

-0.24 

0,885 

42,6 

0,956 

52,6 

1705 

-0.06 

0.950 

46,7 

1,215 

59,7 

2505 

0.04 

1,115 

27,1 

2,050 

31,6 

25  July  1956 

0405 

0.05 

1.100 

16,9 

1,525 

12,7 

6  Aur  1956 

2105 

0.07 

1,200 

12,1 

1,210 

11,7 

7  Auf  1955 

2005 

0,05 

1.100 

22,6 

1,517 

U.  -  O  N/ 

27  Aur  1956 

2205 

0.04 

1,00 

10,6 

0,386 

7,93 
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10.  Summary  and  conclusionSo 

A  functional  relationship  betv/een  C/Cg  and  Ri  vjas 
obtained  by  usinf  an  expression  for  Cg  sufpested  by  LettaUo 
Better  agreement  was  obtained  after  making  some  revisions 
in  Lettau’s  formula  for  Cg.  An  expression  for  has  the 
practical  significance  that  one  can  get  directly  the  ratio 
of  a  windspeed  at  some  anemometer  height  to  a  surface  geo- 
strophic  wind  by  knowing  the  roughness  parametero 

By  using  the  revised  expression  for  Cg,  a  new  functional 
relationship  between  C/Cg  and  Rl  was  obtainedo  This  function 
has  different  analytical  forms  in  unstable  and  stable  condi¬ 
tions. 

Another  way  of  obtaining  C/Cg  was  based  upon  computa¬ 
tion  of  U-Ht  using  Blackadar’s  results.  By  this  means  a 
v;ind-profile  expression  was  obtained  which  varies  v;ith 
stability.  Prom  this  new  procedure  Dne  can  ret  the  value 
of  Cf/Cg  without  using  C/Cg, 

rl 
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APPENDIX 


Computation  for  Cg  usinf  Lettau’s  formula 


The  orifinal  equation  is 


o.‘oA 


(Ca,  Re»)  —  a.Sij 

By  cross  multiplying 

LgCCc^Ro)  ^  +  o.vo/4/c*^ 


or 


Cc^Ro  =  »o 


Ci.^A  +0  6o4/Ca) 


=  173. S  X  10 


oao4/Cc 


th  en 


r  D  /  oao  4  X  ^-303  \ 

SfJlo  -  exp  ( - ^ - ) 

»73>5  ^  ' 

r  exp  (  O'i'iy  y/Co.) 

Let  Oe  2395/Ca  :::  y  substi tutinr  in  the  above  eov^ation 

T=CV)  =  -e":.  o 

73.5  \y/ 

__  0>239g/ 173. g  Ro  -  y  _  Q 

y 

Since  V  O 


f(y)  «  o.  2335/173.- ffRo 

Takinr  the  lops  to  the  base  10 

£o3(2i^'Rc)=  L.,y  -1-  (48) 

Let  lof  -  constant,  and  replace  lop  y  by 

0^4Z4zh^y  and  substitute  in  (48) 

o.43X,3^^  +  o-Ay^iy~  K  =  © 

To  pet  the  root  of  ecuation  (49),  we  v/ill  use  the 


(49) 


Newton- Haphson  method 


M  h  =  V.-  - 


(50) 


Since  ffy)--  0. 4S43/y +•  Oo 4243,  therefore  by  substltutinf  this 
in  (50)  we  obtain 

y,v,  =  y  -  -t  -  i< 

Substitxitinr  log-  <^ej  for  K  and  divide  the  fraction 

by  0,4343  , 

0.4342 

Vi.,  =  y,  -  -■^V' 

Chanfing-  to  2,3026  log- 


y.  =  y  - 

•'l  +  i  ■'> 


y,”  +  yt'  +6*53^7  “ 


»/y,-  + » 


(51) 
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